Unsteady Hall Magnetohydrodynamics near a Hyperbolic Magnetic 

Neutral Line: An Exact Solution 



Bhimsen K. Shivamoggi* 
Los Alamos National Laboratory 
Los Alamos, NM 87545 



Abstract- Unsteady Hall Magnetohydrodynamics (MHD) near a hyperbolic magnetic 
neutral line is investigated. An exact analytical solution describing a self-similar evolution is 
given. This solution shows a negligible impact on the current-sheet formation process near 
the hyperbolic magnetic neutral line at small times by the Hall effect but, subsequently, a 
quenching by the Hall effect of the finite-time singularity exhibited in ideal MHD and, hence 
a prevention of the current density blow-up at large times. The asymptotic result given by 
this time-dependent solution is in full quantitative agreement with the formulation of steady 
Hall MHD near a X-type magnetic neutral line (Shivamoggi [23]). The latter formulation 
showed that this asymptotic result indeed corresponds to a hyperbolic configuration of the 
magnetic field lines in the steady case. 

1 Introduction 

When a plasma collapses near the neutral line of the applied magnetic field the continual 
accumulation of the magnetic flux in the region of the neutral sheet puts the current sheet in 
a non-stationary state (Syrovatskii [1]). An exact self-similar solution of the MHD equations 
for a time-dependent, two-dimensional (2D) flow of an incompressible plasma in a hyperbolic 
magnetic field was given by Uberoi [2] and Chapman and Kendall [3]. This solution had 
an initialy current-free magnetic field, so it is not appropriate for the reconnection problem. 
This solution was modified (Shivamoggi [4]) so as to relax this restriction, and hence make 
it suitable for the reconnection problem. This solution was generalized to incorporate a 
uniform shear-strain rate in the plasma flow (Shivamoggi [5], Rollins and Shivamoggi [6]), so 
the magnetic field lines now undergo not only sweeping but also shearing by the plasma flow. 
The above solution predicted a sequence of events associated with the evolution of a current 
sheet in a hyperbolic magnetic field, in agreement with laboratory experiments (Frank [7]) 
on the collapse of a plasma near the hyperbolic neutral line. In recognition of the numerical 
results (Brunnel et al. [11]) showing the significant effect plasma density variations near 
the magnetic neutral point have on the magnetic reconnection processes taking place there 
the above solution was generalized further to incorporate density variations in the plasma 
(Rollins and Shivamoggi [12]). The current-sheet formation process was found to speed up 
in the presence of a plasma density build-up near the current sheet. 
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Fast magnetic reconnection processes in laboratory (ex: sawtooth collapse in tokamak 
disclianges) and space (ex: solar flares and magnetospheric substorms) can be described using 
collisionless plasma models (Yamada et al. [13], and Sliibata [14]). In a high-/? collisionless 
plasma, on length scales shorter than the ion skin depth di, the electrons decouple from the 
ions and the electron dynamics is governed by Hall currents (Sonnerup [15]). The decoupling 
of ions and electrons in a narrow region around the magnetic neutral point allows for rapid 
electron flows in the ion dissipation region and hence a faster magnetic reconnection process 
(Mandt et al. [16] and Biskamp et al. [17]). 

It may be mentioned that even in the absence of the Hall currents fast collisionless mag- 
netic reconnection has been shown to be feasible (Scholer et al. [18], Jaroschek et al. [19], 
Bessho and Bhattacharjee [20], Ishizawa and Horiuchi [21]). This is caused by - 

• inductive electric flelds generated by lower hybrid drift instability, or 

• non-gyrotropic pressure tensor effects caused by the ion-meandering motion near the 
magnetic neutral point. 

In recognition of the important role played by the Hall effect in fast magnetic reconnection 
processes an investigation of unsteady Hall MHD near a hyperbolic magnetic neutral lin(S 
is therefore in order - this is the objective of this paper. The asymptotic result given by the 
time-dependent solution in question turns out to be in full quantitative agreement with the 
formulation of steady Hall MHD near a X-type magnetic neutral line (Shivamoggi [23]). 

2 Governing Equations for Hall MHD 

Consider an incompressible, two-fluid, quasi-neutral plasma. The equations governing this 



plasma dynamics are (in usual notation) - 

(9v 1 
^^A-Q^ + (ve ■ V)ve] = -Vpe - ne(E + -Ve X B) + ner^J (1) 

X 

nmi[—^ + (vi • V)vi] = -Vpi + ne(E + -v^ x B) - ne?7J (2) 
ot c 

V • Ve = (3) 

V • Vi = (4) 
V-B = (5) 

V X B = -J (6) 



^This problem has also recently been considered by Litvinenko [22]; however, his solution 
is rather different from the present solution and does not show the asymptotic connection 
with steady Hall MHD near the X-type magnetic neutral line. 
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where, 



J = ne(vj - Ve). (8) 

Neglecting electron inertia (mg =^ 0), equations (1) and (2) can be combined to give an 
ion-fluid equation of motion - 

1 

^"^^[^ + • V)^*] = -^(Pe + P») + - J X B (9) 

ot c 
and a generalized Ohm's law - 

1 1 

E + -ViXB = r/J + J X B. (10) 

c nec 

Non-dimensionalize distance with respect to a typical length scale a, magnetic field with 
respect to a typical magnetic field strength Bq, time with respect to the reference Alfven 
time ta = where V^o = ^ ' ^ = rriin, and introduce the magnetic stream function 
according to 

B^Vi/j xU + bl, (11) 

and write the ion-fluid velocity as 

Vi = (i^ X V,) xi^ + wU = v + wi^ (12) 

and assume the physical quantities of interest have no variation along the 2;— direction. 
Equations (9) and (10) then yield 

+ (v • V)]v = -VP - (V^V) VV' (13) 

[^ + (vV)]V' + a[6,V']=77V'V' (14) 
+ (v ■ V)]6 + a[^, W'^] + [i/j, w] = (15) 
[^ + (v-V)]^/;-[6,^] = (16) 



where, 



a 
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3 Hall MHD Near a Hyperbolic Magnetic Neutral Line 



Consider the initial- value problem near a hyperbolic magnetic neutral line in Hall MHD with 
initial conditions - 

t = : = -7oX, Vy = joy, w ^ -1/ a{kx^ - y^), ip ^ kx'^ - y'^ , b^Cxy (17) 

where 70 and k are externally determined parameters with 70 > and C > 0. This initial 
condition describes a stagnation-point plasma flow impinging transversely onto the a; = 
plane and incorporates equations (4) and (5). The spatial structure for the out-of- plane 
magnetic field described by this initial condition is in recognition of the quadrupolar out-of- 
plane magnetic field h pattern characterizing the Hall effects (Terasawa [24]). Laboratory 
experiments (Ren et al. [25]) and in situ measurements in the magnetotail (Fujimoto et 
al. [26], Nagai et al. [27], Oieroset et al. [28]) have also confirmed the latter signature of 
the Hall effect. The Hall magnetic field h is believed to be produced by the dragging of 
the in-plane magnetic field in the out-of-plane direction by the electrons near the X-type 
magnetic neutral line ([16]). 

The Lorentz force due to the initial magnetic field is 

i = : J X B = 4(1 - k)kxi^ - 4(1 - k)y\y - 2C{kx^ + y^%. (18) 

Wc take k > 1 so that this Lorentz force is directed so as to maintain the prescribed initial 

stagnation-point fiow. 

Let us assume that the solution, for t > 0, of equations (4), (5) and (13) - (16) with the 
above initial conditions is of the self-similar form 

Vx(x,y,t) = -^(t)x, Vy{x,y,t) ^ ^{t)y 

w{x,y,t) = -[f3{t)y^ - ka{t)x^] 
a 

ijix,y,t) = ka{t)x^ - (3{t)y^ 
V: t) = Cxy 

P{x,y,t) = -^u{t){x'' + y^) + Po, iy{t) > (19) 

with 

t = = = 1,7 = 70. (20) 

For the solution (19), V^V' = f{t) and — 0, so the effect of resistivity in this case is to 
add a function of t to (which leaves the magnetic field unaltered) and hence to introduce 
an electric field along the z-axis. Wc therefore drop the resistivity in the following. Further, 
for an incompressible plasma, the pressure does not have a dynamical role. So, it is forced 
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to be an enslaved variable in the sense that its form is chosen so as to be compatible with 
equations (13) - (16) given the ansaetze for Vx,Vy,w,ilj and feUH 
Substituting (19) into equation(13), we obtain 

7 = 2(A;V-/?2) (21) 

while equation (14) gives 

d - 2(7 + aC)a = (22) 

/3 + 2(7 + aC)f3 = (23) 

Equations (15) and (16) are identically satisfied by the solution (19). 
We have from equations (21)-(23), 

a{t) = e2(^+-^*) (24) 

P(t) = e-2(^+-^*) (25) 
^^ + 2aC^= [/tV(^+"^*) + e-^(^+"^*)] - A (26) 

with 

t = : 7 = 0. (27) 
Equation (26) along with (20) and (27) yields 

A = {k^ + 1)- 72 - 2aC7o. (28) 

It may be noted that (24) and (25) are consistent with the ion-fluid incompressibility condi- 
tion - 

^ (-in [a{t)P{t)]] = V • V = (29) 



dt \2 



2ln the generic situation, it may be mentioned that there are six parameters - two for 
the magnetic flux function ip, one for the out-of-plane magnetic field b, and three for the 
velocity field (v, w), to be determined by only five scalar equations (13)-(16), so there is 
non-uniqueness in the solution. This is resolved by specifying the parameters as in (19) to 
close the system, so (19) is one exact solution. However, this exact solution turns out to 
have considerable physical significance, as discussed in the following. 

3lt may be noted that Litvinenko [22] sets up the solution and the inital conditions for 
w and b rather different from those prescribed in (19) and (20). In Litvinenko's solution 
the amplitude of the driving plasma flow velocity field is taken to be time indepentent, 
i.e., 7 = const, while in the present solution the out-of-plane magnetic field is taken to be 
time independent, i.e., b = Cxy,C = const. This difference turns out to lend considerable 
physical significance to the present solution. 



5 



which is derivable from equation (14) on substituting (19). 
For small t, equation (26) gives 

7(i) ^ -fot + {k^ - l)t^ (30) 

while for large t, equation (26) gives 

7(t) —aCt or 'y{t) ^ —aC, t large (31 a,b) 

(30) shows that the Hall effect (a ^ 0) does not materiahze to 0{t^). However, (31) 
shows that the finite-time singularity exhibited in ideal MHD (Shivamoggi [29]) is quenched 
by the Hall effect. Thus, though the Hall effect does not impact the current-sheet formation 
process for small times, it prevents the current density blow-up at large times. This result 
may be further appreciated by noting that equations (24)- (26) lead to the exact invariant - 

(7 + aC f - {k^a^ + (3^) = const. (32) 

(32) clearly shows the suppression of the plasma collapse process near a hyperbolic magnetic 
neutral line in Hall MHD, for large times (when (31b) becomes valid). Physically, the 
suppression of the plasma collapse process near a hyperbolic magnetic neutral point in Hall 
MHD appears to be caused by the dispersive activity of whistler waves which is known to lead 
to current-sheet broadening, as confirmed by laboratory experiments (Urrutia ct al. [30]) 
and satellite observations at the magnetopause and the magnctotail plasma sheet (Sonnerup 
et al. [31], Fairfield et al. [32]) as well as numerical simulations (Shay et al. [33] and [34]). 
(31b), in conjunction with (19), also shows that, for large t, the level curves of the out-of- 
plane magnetic field are also the streamlines of the in-plane ion flow. It is pertinent to note 
that the asymptotic result (31b) is in full quantitative agreement with the formulation of 
steady Hall MHD near a X-type magnetic neutral line ([23]). The latter formulation showed 
that (31b) indeed corresponds to a hyperbolic configuration of the magnetic field lines in the 
steady state. 

4 Discussion 

In recognition of the important role played by the Hall effect in fast magnetic reconnection 
processes this paper makes an investigation of unsteady Hall MHD near a hyperbolic mag- 
netic neutral line in Hall MHD. The Hall effect is found not to impact the current-sheet 
evolution process near the hyperbolic magnetic neutral line for small times. However, sub- 
sequently, the Hall effect is found to quench the finite-time singularity exhibited in ideal 
MHD and hence to prevent the current-density blow-up at large times. The asymptotic re- 
sult (31b) is in full quantitative agreement with the recent formulation of steady Hall MHD 
near a X-type magnetic neutral line ([23]) which showed that (31b) indeed corresponds to a 
hyperbolic configuration of the magnetic field lines in the steady case. Besides, in this range 
of time, the level curves of the out-of-plane magnetic field are also the streamlines of the 
in-plane ion flow. 
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